Hadamard Re normalization of the Stress Energy Tensor on the Horizons of a 
Spherically Symmetric Black Hole Space-Time 
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We consider a quantum field which is in a Hartle-Hawking state propagating in a general spher- 
ically symmetric black hole space-time. We make use of uniform approximations to the radial 
equation to calculate the components of the stress tensor, renormalized using the Hadamard form 
of the Green's function, on the horizons of this space-time. We then specialize these results to the 
case of the 'lukewarm' Reissner-Nordstrom-de Sitter black hole and derive some conditions on the 
stress tensor for the regularity of the Hartle-Hawking state. 
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I. INTRODUCTION 

The renormalized expectation value of the stress en- 
ergy tensor operator {T^i,)ren is of fundamental impor- 
tance in the study of quantum field theory in curved 
space time as it governs, via the semi-classical Einstein 
field equations 



I ren i 



(1.1) 



the back-reaction of the quantum field on the geometry 
of the space-time. There is a long and fruitful history of 
calculations of {Tfj_^)ren for various space-times, begin- 
ning with the work of Candelas and Howard for a mass- 
less, conformally coupled scalar field in the Schwarzschild 
black hole space-time [H-Q. This was then extended to 
the case of scalar fields with arbitrary mass and cou- 
pling to the Ricci scalar, in general spherically symmet- 
ric space-times, by Anderson, Hiscock and Samuel 0, 
Other examples of calculations of {T^u)ren can be found 
in [^13. 

These calculations relied on the renormalization coun- 
terterms, denoted by {Tu,i,)ds^ which were first calcu- 
lated by Christensen [13[ . This calculation in turn relied 
on the DeWitt series representation for the Green's func- 
tion, which is an asymptotic power series in inverse pow- 
ers of the mass of the field m. This series is ill-defined for 
the massless case and requires some severe modification 
in order to be applicable to this case. We choose instead 
to follow the Hadamard renormalization procedure in the 
formulation of Brown and Ottewill jT3 | which, building 
upon the axiomatic approach of Wald jTH, constructs 
a renormalization prescription which is well defined for 
both massive and massless fields. 

The calculation of the components of (T^jy)ren in the 
exterior region of a black hole space-time naturally splits 
into two distinct parts, a quasi-analytical calculation to 
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find the exact values on the horizons and a numerical 
calculation which is valid in the exterior region excluding 
the immediate vicinity of the horizon. This paper is con- 
cerned with the former calculation, while in an upcoming 
paper (henceforth referred to as Paper H) , we will discuss 
how the application of the Hadamard renormalization 
procedure leads to an alternate method to that of An- 
derson, Hiscock and Samuel f^, for calculating {Tfii,)ren 
in a general spherically symmetric space-time. 

We note that these calculations do not assume that the 
exterior region of the space-time possess a single hori- 
zon, in fact, once the general calculations are done, we 
will apply the results to the case of a lukewarm black 
hole, whose exterior region possesses both an event and 
a cosmological horizon. 

In a previous paper , we developed a technique for 
calculating the vacuum polarization of a scalar field on 
the horizon of a spherically symmetric black hole space- 
time. This paper was motivated by the numerical can- 
cellation of horizon divergences to give a finite answer, 
which arose in a paper by Winstanley and Young on the 
vacuum polarization for lukewarm black holes [17| . It is 
therefore natural to ask if we may extend the technique 
developed in to the stress tensor case, allowing us 
to obtain finite values on the horizon without relying on 
numerical cancellations of divergent quantities. 

This paper is organized as follows, in Sec. |ll] we will 
outline the Hadamard renormalization procedure. In 
Sec. mil we will demonstrate the application of this 
method to obtain unrenormalized mode sum expressions 
and their renormalization subtraction counterterms. In 
Sec. IIVI we will calculate these mode sums using a uni- 
form approximation to the radial solutions allowing us to 
obtain renormalized values which we then apply to the 
lukewarm case in Sec.|Vl Finally our conclusions are pre- 
sented in Sec. |Vll Throughout this paper we use the sign 
convention of Misner, Thorne and Wheeler [l^ and we 
will work in units in which SttG = h = c = = ^- 
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II. HADAMARD RENORMALIZATION 

In this paper we will be concerned with the calculation 
of {T^jj) resulting from quantum scalar fields. The action 
for a scalar field (f) with mass to on a background space- 
time of dimension n, possessing metric g^n^ is given by 



S 



-9 



(Px{g^'''cp,^(j),, + [m2 + ^i?]02}^ (2.1) 



where g denotes the determinant of g^^, R is the Ricci 
scalar and ^ is a dimensionless constant describing the 
coupling between the scalar field and the gravitational 
field. Requiring that the variation of the action with 
respect to (j) vanishes yields the scalar field equation 



(□- 



(2.2) 



where □ = g^"V ^^Vu- The classical stress tensor is de- 
fined by the equation: 



(1 - 2O</'^^0^'' + (2e - ^)5^''</';a</'^^ 

+ 2£,g^"'(t)U(t> -f S,{R^"' - \Rg^''')(t>'^ 



(2.3) 



Vacuum expectation values of the various quadratic 
products of field operators can be identified with var- 
ious Green's functions of the wave equation (|2.2p . Of 
particular importance in the study of quantum field the- 
ory in curved space-time is the Feynman propagator Gf, 
defined as the time ordered product of the fields [l^ 

Gf{x,x') = t{0\T{${x)${x'))\0), 

=i{e{t - t'){Q\Ux)^{x')\Q) + e{t - t'){Q\^{x')^{x)\Q)). 

(2.4) 

The Feynman propagator satisfies the equation 

5"(x,x') 



{U-m^ -^R)Gf{x,x') 
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(2.5) 



where (5"'(x,x') is the n-dimensional delta function. 

We will perform our calculations on a static spheri- 
cally symmetric black hole background spacetime with 
line element: 



ds2 ^ -f{r)df 



1 



fir) 



■ si-n^ 9d(j)'^ . (2.6) 



The spacetime will have a horizon at r = tq whenever 
fi'To) = and in this case the surface gravity of that 
horizon is given by kq = |/'(ro)|/2 . Following the stan- 
dard procedure we Euclideanize our spacetime, that is we 
perform a Wick rotation t — ^ it, Eq. (|2.6p then becomes 



ds2 = f{r)dT^ + 



1 



fir) 



dr^ + r^de^ 



2 • 2 

r sm 



ed(f^. (2.7) 



Assuming kq 7^ 0, this space will have a conical singu- 
larity whenever /(rg) = which may be removed by 
making r periodic with period lixjnQ. This periodicity 
in the Euclidean section corresponds in quantum field 
theory on the Lorentzian section to a thermal state with 
temperature T — kqII-k. 

In this configuration the quantum field satisfies a 
wave equation where □ is now the Laplacian on Euclidean 
space. The Euclidean Green's function Gsix, x'), for this 
elliptic equation, is related to Gf by 

Gpit, x; t\ x') = -iGsiir, x; ir', x'), (2.8) 

and satisfies the equation 

{a-m'-^R)GEix,x') = -^-^^^, (2.9) 

with X = {T,r,d,ip). A key advantage of Euclidean field 

theory is that □ is now an elliptic operator and hence 

has a unique, well defined inverse when supplemented by 

appropriate boundary conditions [20j . 

We may now define a formal expression for (T^i/) given 

by: 



(f^"^) =n[\imT^"'GEix,x')], 



(2.10) 



where TZ denotes taking the real part, and r'^'^ is a dif- 
ferential operator which reduces to Eq. (|2.3[) in the coin- 
cidence limit, for example (l^ 

r'^"^ = (1 - 20.9.rV^V''' + (2^ - i).9^^g„?V"'V„ 

-2fV^V'' + 2eg^'^V„V" + ^{R'"' - ^Rg"") - ^m^g"", 

(2.11) 

where gjf is the bivector of parallel transport, which 
serves to parallel transport a vector at x' to a vector 
at X. 

These expressions are perfectly valid in the classi- 
cal theory, however in the quantum theory, being con- 
structed from products of distributions evaluated at the 
same point, they are divergent. The characterization of 
these divergences (and their identification as an adjust- 
ment to the parameters of the theory) is the role of renor- 
malization theory. As mentioned in the introduction the 
usual approach to renormalization begins with the De- 
Witt series representation. We choose to adopt an ap- 
proach based upon the Hadamard series representation 
to the Euclidean Green's function, as it is equivalent to 
the standard method but is better defined, in particular, 
for the massless theory [l3]. We will now briefly outline 
this method. 

In four dimensions the Green's function possesses 
the following singularity structure, first identified by 
Hadamard ^21j, ^2^] : 



GEix,x') 



1 

8^ 



Ai/2 



+ Vln{Xa) + W 



(2.12) 
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where 2a(x,x') denotes the square of the geodesic dis- where vi = [Vi{x,x')] and is given by the expression [1^ 



tance between x and x' , and 
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(2.13) 



is the biscalar form of the VanVleck-Morette determi- 
nant. V{x, x') and W(x, x') are regular biscalar functions 
possessing expansions of the form 

oo oo 

V{x, x') = Mx, a;')CT"; Wix, x') = ^ W„(x, x')fT", 



(2.14) 

where Vn{x, x') and Wn{x, x') are themselves regular bis- 
calar functions. A is a constant which is required to en- 
sure that the argument of the log function is dimension- 
less. Following the convention of Christensen [l^, we let 
A — where /i = m for a massive field but is ar- 

bitrary when the field is massless Following [5], we 
choose to set fj, = 1 for massless fields. 
Imposing Eq. (|2.9p for x x' , and demanding balance of 
explicit powers of a. one can derive recursion relations 
for the coefficients Vn{x,x') for n > and for Wn{x,x') 
[2^ for all n > 0. Wo{x,x') remains undetermined and 
corresponds to the freedom to add to Ge{x, x') solutions 
of the homogeneous wave equation. The Vn{x,x') are 
purely geometrical, all the information about the state 
in contained in W{x,x'). All the singular behaviour of 
the Green's function is geometrical and state indepen- 
dent. 

We return now to our definition of the stress tensor 



{f^"') = lim t^"'Ge{x,x'), 



(2.15) 



where r^" is defined in Eq. (|2.1ip and it is implied that 
we are only taking the real part of the limit. As we have 
previously stated, this expression is divergent and is a 
priori meaningless. However, as the singular behaviour 
is geometrical, the difference between the stress tensor 
of two different quantum states, \A) and \B) say, is well 
defined and given by 

{A\f'"'\A) - {B\ft"'\B) = 
lim {t^'''Gea{x, x') - rf^'GEBix, x')} 

x—^x' 

= [tI"'Wa{x, x') - t'"'Wb{x, x')l (2.16) 

where the square brackets denote that the coincidence 
limit has been taken. Now the tensor 



r^^[W^^] = K^W^A(x,x')] 



(2.17) 



is well defined and contains all the state dependent in- 
formation. We will use it to construct a new definition 
for the stress tensor. Firstly we remark that although 
t'^'^[Wa\ is symmetric it is not, in general, conserved. In 
fact it can be shown that (l3 [ 



-2v'f, 



(2.18) 



+ \i^-WR' + \rr^'i^-\)R+'^- 

(2.19) 

If we therefore choose to define our renormalized stress 
tensor for a quantum state \A) by 



{fnren = ^ {t^'" [Wa] + ^VlQ^ , (2.20) 



then we see that our definition of the stress tensor is 
conserved and has trace [l3| 



(V) = i (2^1 + Wi- mwA{x) - m^WA{x)) . 

(2.21) 



where wa = [W^A(a:, a:')]. In the conformal case this defi- 
nition reproduces the standard trace anomaly vi/An'^. 
Wald has shown that if a given definition for a stress ten- 
sor satisfies certain natural conditions then it is unique 
up to the possible addition of conserved geometrical ten- 
sors corresponding to the metric variation of the terms in 
the gravational action [l^ . The definition (|2.20p satisfies 
these axioms 



14{ and so is equivalent to any (T^,y) 



derived using one of the other renormalizing techniques, 
but has the added advantage of being more direct. As 
noted above we must introduce an implicit length scale 
in order to make the argument of the logarithm func- 
tion in Eq. (|2.12p dimensionless. A different choice in 
length scale results in a different definition for {Tn^)ren, 
however these tensors will differ only by a multiple of 
[r'-'-'^V{x, x')], which is a conserved tensor in Wald's class. 

Mc Laughlin has shown that renormalization using the 
Hadamard regularization procedure is equivalent to co- 
variant geodesic point separation technique In fact 
he has proven that, for scalar fields, the stress tensor ob- 
tained via Christensen's covariant geodesic point separa- 
tion technique procedure is equal to that obtained using 
the Hadamard form plus a geometric term 



M 



167r2 



2 U!/ 

-TO 



(2.22) 



which can be absorbed into the left hand side of the semi 
classical Einstein equation in the usual way by renormal- 
izing the constants G and A. ^A^^'^ is clearly a conserved 
quantity, therefore we are free to add it to our definition 
of the stress tensor, without violating Wald's axioms, in 
order for our results to be in agreement with those ob- 
tained using Christensen's method. Hence we choose to 
give a final definition of the renormalized stress tensor as 

{fnren = ^ {t^^^IWa] + ^v^g^ + (2.23) 
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with trace 



(2.24) 



again giving the correct trace anomaly. 
We note here that this analysis has been extended to arbi- 
trary spacetime dimensions by Decanini and Folacci ^23] . 



III. APPLICATION OF HADAMARD 
RENORMALIZATION 

In practice, we calculate r^"^ [Wa] on a given space-time 
by subtracting off the action of r'^'^ on the singular part 
of the Hadamard form (j2.12p from an expression for the 
unrenormalized value oir'^'^GEix, x'), then finally taking 
the coincidence limit. The unrenormalized expression we 
make use of was first derived for a Schwarzschild space- 
time by Candelas [l] and extended to the case of general 
spherically symmetric space-time by Anderson [3] and is 
given by 



Ge{x,x')^£-^ e»«(---')^FKcos7)Xn;(r-,r'). 

(3.1) 



Where n labels the mode's frequency, I labels the mode's 
total angular momentum, and XnK^i^') the Green's 
function for the radial equation 



dr 



dr 



n^nl 1(1 + I) n , , 



5{r- 



(3.2) 



We define pni and g„/ to be the independent solutions to 
the homogeneous version of this equation, with pni de- 
fined to be the solution which is regular on the lower limit 
of the region under consideration, while qni is regular at 
the upper limit of the region. Xn/(r, r') is then given by 



Xni{r,r') = CniPni{r<)qni{r>), 



(3.3) 



where r< is the lesser of the two points r, r' , and r> the 
greater. Cni is fixed by the Wronskian condition 



Pni- 



dr 



qni- 



dpnl 



dr 



1 



127- (3.4) 



Since we have that Xni{r, r') is real, we may then express 
Ge{x,x') for a thermal state in a form which we will use 
for the remainder of this paper: 



Ge{x,x) 



n=0 



F{n) cos(n/v(e^))^(2/-t- l)Pi(cos7) 

1=0 

xC„ip„i(r<)g„;(r>), 
(3.5) 



where = t - r', F(0) = k/Stt^ and F{n) = kI^-k^, 
n > 0. Henceforth we choose to set r> = r. In static, 
spherically symmetric space-times, states respecting the 
same symmetries have a stress tenor {T'^^)ren, which is 
diagonal. Using the definition given in the previous sec- 
tion we may express the diagonal elements of (r'^)j.en in 
the following manner: 

{F jy} ren "^(2 G -lyi/^^ren + (2^ 2)!-^ G •act']re7 



2e[5""G: 



aa] ren 

(3.6) 



where v is not summed over, ; denotes covariant differ- 
entiation and we have dropped the {x,x') dependence as 
well as the E subscript for notational convenience. 

Our strategy for computing these components of 
{T^^)ren IS as foUows. Firstly we calculate the required 
derivatives of Ge{x, x'), after which we may take the par- 
tial coincidence limit {t — > t' ,6 — > 6*', </) — > </>'} and also we 
place r' on r^ . We then rcnormalizc these expressions by 
subtracting off their respective singular Hadamard parts 
before taking the coincidence limit r — > rg. Finally we 
insert these renormalized expressions into Eq. p.6|) . 



A. Unrenormalized mode sum expressions 

Before we can begin to consider working with Eq. 
we need to first calculate these bivectors for our partic- 
ular choice of point separation. The bivector of parallel 
transport are defined by the equation 



O" gab' -.a' 



0, 



(3.7) 



with the boundary conditions that gab' — gab when 
X — x' . For a general point separation this expression 
can be quite complicated, however for radial separation 
this reduces to quite a simple form. By virtue of the 
symmetries of the space-time we have 



hence Eq. p.7p reduces to 

CT''' gab';r' — 0. 



(3.8) 



(3.9) 



Using the definition of the covariant derivative we see 
that the bivectors of parallel transport for radial separa- 
tion are determined by 



gab'.r' — b'r'go-p' 



(3.10) 
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We require the non- vanishing components of the connec- 
tion, these are 



f'jr) 
VirV 



Tle = -rf{r); T^. = -r/(r) sin^ 



r;* = ^^^r«, = C = i; r^, = -sin0cos( 



r<l> — 4>r — ^ 90 — 09 — ^ tr — 



_ fir) 
2/(r)- 
(3.11) 



The calculation of the bivectors is repetitive in nature so 
we will just show the calculation of the one and the rest 
follow similarly. 
For Qrr' we have: 



9rr' ,r' — ^ r'r' 3rp' — ^r'r'dr 



fir') 



2/(r') 



■j-grr'- (3.12) 



Integration gives 



fir) 



1/2 



Jir')J VfWF) 
The other components are given by: 

'firV'^' 



(3.13) 



gw 



fir) 



9tt 



= -V/W/(r'). (3.14) 



966' = —gee = rr 
r 



= rr' sin^ ( 



(3.15) 



(3.16) 



Where we have used the relation d/dt = id /dr. Taking 
the partial coincidence limit yields 

oo oo 

{9''G.,u} = E E(2^ + ^)Pniir')qniir) 

(3.19) 

This vanishes in the limit r ro as Pniiro) = for n > 

as 

Pniir') = ao(r' - ro)"/^ + O ((r' - ro)"/2+i) , (3.20) 
with ao = 1/ \//tro • Likewise taking r' — )• ro 



1=0 



dr 



(3.21) 



and hence 



fir) apK 
2 47r2 



(3.22) 



Repeating this procedure for the other derivatives leads 
to the following mode sum expressions 

oo 

{Ge} = ^ao ^i2l + l)qoiir). (3.23) 



;=o 



{g^'G.,,} ' ^ "° 



^(2Z + l)gu(r) (3.24) 



It is straightforward to show that all other components 
of Qab' vanish. 

We turn our attention to calculating the derivatives of 
Gsix, x'). For the purposes of this calculation we choose 
to denote the partial coincidence limit, {t t',9 — 
6', (j) — >■ (/>'}, of a bitensor A{x, x') by {A}. We will show 
the derivation details for one derivative {G;tt}, the oth- 
ers follow along similar lines. Using the definition of the 
covariant derivative, we have that 

G;tt = G^tt - TuG,,. (3.17) 

We consider the partial derivative first 

92 ^ 



G,« = -^:^Eco«( 

X Pniir')qniir) 



nK€r )2^(2; + l)P((cos7) 

1=0 



= -^-^ ri'r? cos(nKex) ^^(2/ + l)P;(cos7) 

n=0 1=0 

xp„j(r')g„i(r). (3.18) 



{.-'G..} = :^.AA^aoV(2l + l)^ (3.25) 



4^2 V 2 ""^^ 

(=0 



{G-es'} _ {G;00'} _ Koo x^icy. 

r2 ~r2sin2^~47r2£^^^' + '^V 



dr 



fir) dqoi{r) 



r dr 

i{i + r 



2r2 



-qoiir) 



(3.26) 



{G-66'} _ {G;00'} i^ap ^ + 

^~2l{ = 1 2 / ^v2t + 1) — 7: Qoi[r) 

rro rrosm^e 4ro7r2 f-' ' 2r 



(3.27) 



Note that we did not calculate a mode sum expression for 
fir){G-rr} since, as will be seen shortly, we may exploit 

the wave equation to obtain [g^''Grr]ren once we have cal- 
culated the renormaized value of the other components. 



6 



B. Renormalization subtraction terms 

As discussed earlier, the renormalization subtraction 
terms for a particular component are obtained by the 
action of the corresponding differential operator on the 
singular part of the Hadamard form of Ge, namely 



1 



+ yiog(Aa) 



(3.28) 



Using the expansion 



cr = iga/jAx^Aa;'^' + ^0^3^ Ax" Ax'' Ax'^ 



"Ax'^ +^a;3^Ax"Ax'^' 
+B„^^5Ax"Ax'^Ax'^Ax'^' + . . . (3.29) 



with 



1 



(3.30) 



Babcd ^ ^ 3 ( ^(abc.d) + ( -^9{ab,\a\A\l3\cd) 



2^a(afc^|^|cd) 



(3.31) 



(the higher order terms are easily calculated by hand or 
by using Mathematica Wardeh and Ottewill |26|] 

have developed a Mathematica notebook which allows 
one to expand both A^/^(x,x') and ^(x, x') in a coor- 
dinate expansion in powers of the coordinate separation 
of Ax", for arbitrary point splitting to high order. This 
notebook thus allows one to obtain a series expansion of 
Gsmg(x,x') in powers of Ax", up to the required order 
to capture both the divergence and the finite remainder 
terms of Gsingix, x') and its derivatives in the coincidence 
limit. 

This method works perfectly for any regular point of the 
metric, indeed we will make use of this method in paper 
II to calculate the subtraction terms for temporal sep- 
aration. In the case under consideration here, namely 
the case of radial separation about the horizon, we must 
modify this method slightly. This is due to the singular- 
ity, at r = ro, of this coordinate system. This singularity 



is manifested through the function f{r) which vanishes 
on the horizon, hence Eq. p.29p . which is an expansion 
of a for small e = r — and fixed /(r), is no longer valid. 
To overcome this, we use the definition of a in terms of 
proper distance s: 



2a 



(3.32) 



for space- like geodesies. A great advantage of radial point 
splitting is that it allows one to integrate the line element 
to get an expression for s. Since t = t' , 6 = 6' and (j) = (j)' , 
the line element becomes 



dsi = 



where Sr denotes the proper distance along a radial 
geodesic. Hence 



-.dr'. 



(3.33) 



For a Ricci-fiat space-time (/(r) quadratic), one can per- 
form this integral exactly and hence obtain an expression 
for a for radial splitting everywhere, without recourse to 
the expansion method. Unfortunately, for a general non 
Ricci-flat space-time, this is not the case. We can, how- 
ever, expand the integrand about the horizon, then inte- 
gration yields an expression for in terms of e = r — tq . 
Using this relation we may obtain an expression for a and 
hence A^/^(x, x') and V{x^ x') which are now valid for ra- 
dial separation about the horizon. Implementation of the 
method of Wardell and Ottewill, adapted in the manner 
just outlined, leads to the expressions listed below for the 
required subtraction terms. We note here that for sim- 
plicity, we choose to perform our calculations for a space- 
time with a constant Ricci scalar R. Therfore we may 
define an effective field mass rh = ^Jm? -I- (^ — l/6)i?. 
It is worth noting also that space-times with a cosmo- 
logical constant that are solutions of Einstein's vacuum 
equations, possess a constant Ricci scalar. 



J 



{G si7ig\ 



rrr / Ae 



487r2ro 



0(eln(e)) 



(3.34) 



{U^^G Esing;tt} 



87r2e2 167r2e 



115207r2r4 



-4 (/r + 60/^'m2 - 180m^) + 4} In (e) + 0(e ln(e)), (3.35) 



{g^^GMng-M = {g'^'^G.ing;;^^} = + {R + 6m^) + Fee 

+ ^^g^Q^2^4 {4/0^2 - 2f^rl (2/^' + /^"ro + GOm^) + (/q'" + ISOm^) + I20m\l - 4} ln(e) + 0(e ln(e)) (3.36) 
+ ^^g^Q^,^4 {^o/o (io' + eOm^) + 4roV;' - ^rlf'^ {rofmo ' SOm^) - 4 (45mV^ + l) } In (e) + 0(e ln(e)) (3.37) 

2 -^2 
r tt' /~i 1 /^TTl 



+ 



115207r2r^ 



{^o/o (/o + 60m2) + 4r2/o' -2ri^/o (ro/o")2/o " eOm^) - 4 (45m + l) } In (e) + 0(e ln(e)) (3.38) 



^2 



{-4r^ (/o' + 30m2) - 4 {-V'ofo + f'o + 180m^) + 4/^^o (/o + SOm^) + 4} In (e) + 0(eln(e)), 

(3.39) 
where 

1 r. /A 



^" = 11520^2,4 {in (;^) [-4/f + V'orl {V^ + /^r-Q - eOm^) - ((/^'^ + SOj^'m^ - ISOm^)) + 4 

-4/rr^(ln(2) - 15) + f^rl (4/^'(ln(2) - 3) + 2,Cro(ln(2) - 92) - 120^2(1 + ln(2))) 

+^0 (-/o''(ln(2) - 2) - 60/^'m2(ln(2) - 4) + 180771^(1 + ln(2))) + 4 + ln(16)} . (3.40) 



Fee = {in (A 



4rl + 30m2) + ,-4 (-2/^/^" + /^'^ + 180m^) - 4/^^ (/^' + 30™^) - 4 

+ro [4/^'Vo(39 + ln(2)) - (r^ (/^'(248 + ln(16) + .Cro(4 + ln(4)) + 120m2(ln(2) - 1)) + 280) 

+/o V^(2 + ln(2)) + 60mVo (mVg(3 + ln(8)) + ln(4))] + 4(1 - ln(2))} (3.41) 



115207rM 



o/o (/o + eOm^) + Arlfo" - 2r3/^ (ro/^' + 2/^ - SOm^) _ 4 (45mV^ + l)] 

+4/f r2(2 + ln(2)) + /^'V^ ln(2) + 60/^mV^(4 + ln(2)) - 4(3 + ln(2)) (45m + l) 

-2fo4 [/o (28 + ln(4) + i^'Vo(45 + ln(2)) - 60771^(3 + ln(2))] } (3.42) 

i i^'' {id ^ ^^"^'^ ^ " ^ ~ ^^"^'^ " ^ ^ 

ln(16) + /;'V^(ln(2) - 2) + 60/^mV^(2 + ln(2)) - 4(1 + ln(2)) (45m + l) 

-2/^r^ [/^'(24 + ln(4)) + /^Vo(13 + ln(2)) - 60771^(1 + ln(2))] } (3.43) 



- 1152^ {f^ 



-4rl + 3077.2) _ ^4 (_2/^/- + /„2 ^ 180^4) ^ 4_^/^3 ^ 30^2^ ^ 4 

+ro [-4/^Vo(ln(2) - 21) + {r^ [/ain(16) - 52) + /^'Vo(4 + ln(4)) + 120771^(2 + ln(2))] + 160} 

-^0 (/?r^(2 + ln(2)) + 60771^ (mV^(3 + ln(8)) + ln(4)))] - 4 + ln(16)} (3.44) 
I 

The extension of the above calculation to non-Ricci con- IV. RENORMALIZED HORIZON VALUES 

stant space times is straightforward. 

In this section we introduce uniform approximations 
to the radial function Qni, which will allow us to calcu- 
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late the unrenormalized mode sums to the required order. 
We then expand the resulting expression in the near hori- 
zon limit and show that the divergent terms will cancel 
with those contained in the counterterms calculated in 
the previous section. Finally we will the take the horizon 
limit leaving us with the renormalized horizon values we 
desired. 



A. Approximations of qoi{r) and gi;(r) 

In a previous paper [T^, we demonstrated that the 
EGL uniform approximation to qoi{r) captured enough 
of the horizon behaviour to facilitate the calculation of 
the vacuum polarization on the horizons of a spherically 
symmetric black hole space-time. In order to perform 
the corresponding calculations for the stress tensor, we 
require an uniform approximation which captures more 
of the horizon behaviour. The analysis contained in 
can be extended, giving uniform approximations to both 
Pniif) and qniif) in terms of Whittaker functions. For 
the purposes of this paper, we require only the zeroth 
order approximation to qni{r), which is given by: 



F2{n) 



(2Kr2)«/22i/2^ii-")/4 



n = 



n > 0. 



(4.1) 



Here F is the gamma function, is a Whittaker function 
of the second kind and 




V'« = ^ro[/a-3(r.^-4)r„V^" 
-t-8ro/^ (2n2 - 60^ + 7) + 120 (mVg -f 2^)) 
-8rg/^7o (2«' + - 8) + (4n2 - l) rjl/?] (4.2) 
with Vb = l{l 1) (m2 -f Ci?o)r-o- 

We will now demonstrate that the zeroth order approx- 
imation captures enough of the near horizon behaviour 
of q^i (r) to facilitate the caculation of the required mode 
sums. We will show the details for the n = approxima- 
tion. 

We begin by noting that using standard Frobenius the- 
ory we can find a series expansion for q^i about the event 
horizon. To do this, we require the following expansions 
in e (with Vo/(r) = /(/ + 1) + (m^ + ^R)r^): 



J.Z j (j,) 



Pie + P2e^ + Oie-")] 



(4.3a) 



9ie + g2e^ + 0(e^); 



with 
Po = 1; 
91 



/o 2 

Pi - 7^ + -; 

4k To 



P2 



fll2 

Jo 

16k2 



Vq 



92 = 



,2^2 



^0, 



(4.3b) 



fill 
Jo . 

6k ' 
(4.4) 



The irregular solution obtained by Frobenius analysis on 
the radial equation near the regular singular point r = vq 
then has the form [23 



<loi[r) = -y 



1 - 9ie 



qi{pi 



''^^ - '^V ) ln(., 



(2(?i -Pi)e 
0{e' He)), 



pI-P2- Piqi - Sqf 



_92^2 



(4.5) 



We may form the series for the solution qoi (r) near r — tq 
by adding an appropriate multiple of the regular solution, 
«/Po/ ('-'), to ensure its satisfies the boundary conditions 
at the outer boundary, for example, regularity on the cos- 
mological horizon for a lukewarm black hole or vanishing 
on the outer boundary if the black hole is contained in 
a reflecting box [2^. We note here that since our ap- 
proximation, Q^^^\ is a local approximation, it cannot 
contain all the global information contained in the full 
solution, i.e it will differ form the full solution in its ao; 
term. 

Next, we consider the zeroth order approximation to qoi, 

^W{0) 



'oi 



(r) 



1 



..( 



(er2/)i/4 23/2(^-0)1/4 

It is straightforward to show that (|4.6p satisfies 

d f n „, . d 



2v^oe 



(4.6) 



dr 
with 



[r) 



V^,f(r)QoT°'W=0, (4.7) 



/ 
4 



^2 JI2 



rf") 



+4V'o^(r). 



We now wish to apply Frobenius theory to Eq. (|4.7p . To 
do this analysis we need the equivalent of the expansions 
()4.3ap and (j4.3b|) . Clearly the expansion (j4.3ap is the 
same in this case, so all we need is to find the equivalent 
of Eq. (j4.3bl) for V{r). This is readily computed and we 
find that 



2V^ir) 



r'f{r) 



qie + q2e'^ + 0{e^), 



(4.8) 



with qi and q2 given by Eq. (|4.4p i.e. the potentials Vqi 
and agree to this order. Hence, we can conclude that 
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the local terms of the series expansion of Qolr) about the 
horizon is in agreement with Eq. (|4.5I) up to O [e^ ln(e)) 
and, as discussed above, the global terms ai term will 
differ. So we may then write 

90/ W = QT\r)+l3^Poiir)+n^ir), (4.9) 

where TZ^{r) denotes the remainder terms and is 
0(e^ln(e)) as e — )■ 0. The constants /3qi are determined 
by the requirement that qoi{r) is regular on the outer 
boundary (i.e. contains the correct multiple of the regu- 
lar solution) and are given in Appendix [XI 
We wish to use this approximation to calculate the sum- 
mations 

5i(r)=V(2/ + l)fcllgo/(0, 



S2ir)= fir)} (21 + 1) 



=0 



dr 



(4.10) 



It can shown that for I > e^^/^, the combination 
f^mPoii^) +^oi'('') '^^ts off exponentially in /. Therefore 
we have an asymptotic expression for Si{r) and S2{r) 
valid in the region of the horizon: 



5 



\{r) = £(2/ + l)iii±l)Q„^(«)(.) 0(eln(e)), 



^-0 



oo 



^2(r) = /'(r)^(2? + l)"^«' 

1=0 



dQT\r) 



dr 



0(eln(e)), 

(4.11) 

and so we may conclude that the approximation Q^''^^ (r) 
is of sufficient accuracy for the calculation of the n = 
contribution to the stress tensor. Similar analysis for the 
n ^ 1 approximation Q^^^"^ (r) allows us to conclude that 



quir) = QYi{r)+l3riPuir)+Kir), (4.12) 

where Ru{r) denotes the remainder terms which are 
O (e^/2) as e ^- 0. The constants are determined 
by the requirement that qii{r) is regular on the outer 
boundary and are given in Appendix [X] For the n = 1 
mode we are required to calculate the following sums for 
the stress tensor: 

oo 

V/(0 ,=0 



1=0 



dqiijr) 
dr 



(4.13) 



The argument proceeds along the same lines as in the 
case of go;, allowing us to conclude that 



^3 



^^^(2/ + l)Q]T(r)+0(e), (4.14) 



54 = /ff:(2; + l)^^ + 0(e) 



(4.15) 



1=0 



The calculation of the contributing sums over I of and 
present no problem as they can be shown numerically 
to converge like l^^. 



B. Mode Sum Calculation 

These calculations are tedious and repetitive in nature, 
so for the sake of brevity we will outline the calcula- 
tion details of one component, [g^^ G0g']ren, and simply 
list the results for the_other components. The interested 
reader may refer to [29| for a more detailed discussion. 
We wish to compute the renormalized value, on the black 
hole horizon r = rg, of 



sjs'r \ — ^ 
iff ^-.ee'} — „ 2 2 — 



Stt^Tq r 



(4.16) 



where again, 

S,{r) = £(2/ + l)^l±^Q^(^\r) + 0(eln(e)). 



1=0 



(4.17) 



Unfortunately, we have no way of computing this sum in 
its current form, however we may rexpress it as 



- J2i2l + 1) 



1=0 



1/4 



1=0 

— Sia + Sib, 



1(1 + 1) ( e 



1/4 



K^ikoe'^) 



(4.18) 



where, as above k^ = Vq — ^R^rl + ^. We will now pro- 
ceed to calculate the two component sums in the above 
expression, beginning with as it is the most straight- 
forward. 

a. Evaluation of S'l^ 

We may calculate Sib by making use of the Watson- 
Sommerfeld formula, giving 



Sib = 



r\{\'-\) (^^) ^ Ko{koa\)'/^)dX 



n 



= Ti+l2 



2A(A2 



o27rA 



1/4 



) Kok^{i\)e")d\ 



(4.19) 



Here we have introduced a new integration variable A 
1 + 1/2 and fc^(A) = A^ + m^rl + 1/12 = A^ + kl^ with 
kliiA) 



-A2 -I- P 
A + Kqq. 
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Evaluation of the integral Ii 

We have 

/>oo 

Ii = A{r) / A(A2 - \)KMmrf'^)d\ (4.20) 

JO 

/ , \l/4 

with A{r) = { j . Using the relations 2kQdk = 2XdX 
and kQ{X) — + fcoc '^^ obtain [s^l 



Ii=A{r) X{X' ~ l)Ko{k{X}e/')dX 
Jo 

= A{r) / fco(A)(fc2(A) - fc2^ - \)Ko{ko{X)e/^)dk 

Jko 

(4.21) 



A{r) 



4^2 



Evaluation of the integral I2 

We now consider the integral 



I2 = -n 



2A(A2 + i) / C 



1 + e' 



27rA 



1/4 



i^o(fco(iA)e'/2)dA 
(4.22 



Due to the exponential factor in the denominator, this 
integral is absolutely convergent. Hence we can take the 
near horizon limit inside the sum to give 



I2 = 



1 



-.n 



'0 



2A(A2 + i)dA fl^^^(ko(Ofe ^ 




34 



1 



2 1 1920 V 2 V 2Kr2 



-7^ 



-2A(AM4) 

— — 5^1n(/co(zA))dA 



0(eln(e)) 

(4.23) 



where 7 is Euler's constant and we have used the integral 
formulas [3111 



AdA 



1 + e 



27rA 



1 

48' 



A^dA 



7 



1 + e 



27rA 



1920 



(4.24) 



Note that we truncate the series expansions at the level 
which contributes on the horizon. 

The integral in Eq. (|4.23p has a branch point at A = 
kooiko{iX) — 0), however the divergence is logarithmic 
and therefore is integrable. We then may write 



n 



2X(X^ + -) 

-^—^ln{ko{^X))dX 



'™ ^^^^^^~^\n{ko{^X))dX 



1 + e 
°° 2A(A 



27rA 



2 I 1^ 

4:> 



00 



1 + e 



27rA 



In 



(fco(«A)) 



dX 



(4.25) 



where koiiX) = A^ - /cqq. 
Hence, 

Sib =2^1+2^2 = 

8kl,^K2{kooe/^) - kooe^^Kiikooe^^) 



Air)- 



1 



4^2 

34 / 1 

7 + - In , 

1920 V 2 V 2Kr? 



'^00 2A(A^ + 1) , 

ln(fco(»A))dA 



2A(A2 



1 



o27rA 



In A:o(«A) dA^ +0(eln(e)), (4.26) 



where we haven retained only those terms which con- 
tribute in the e — ?> limit. 

b. Evaluation of the integral Sia 

We calculate Sia in the following manner. Firstly we 
expand the summand in the small e limit and isolate 
the divergence over / by subtracting away the divergent 
terms (if any). We then express the sum over I of these 
divergent terms as a geometric term which diverges in the 
e — > limit and add this back on outside the sum. Fi- 
nally, we apply the Riemann sum argument, developed in 
p^ . to the original summand minus its divergent coun- 
terterm, to give the final value of the sum in the e 
limit. 

Taking the limit inside the sum gives 



Sia — 



1 



{21 + 1) 



[ln(j.o) - + ^0)] 



+0(6) 
(4.27) 

and then, expanding each of the terms in the large / limit, 
we see that 

{21 + l)^^i±ll[ln(.o) - ^(i + -0)] = + 0{l-'). 

(4.28) 

Using the identity derived in [l6j . 

^ 1 



ln(77-l) = -2^ 



(Z + 1)(V^+1)'+ 



- + o(y^). 



and letting rj = r/r^ allows us to express Sia in terms of 
a finite sum over I plus a geometric term which diverges 
logarithmically as e — 7> 0: 



• 1 r (1 + vq) 

{r^f^y/i 23/2(^^0)1/4"'-'^'' 

-A{r)Ko{koe/^) 



W^-.„.o (2v^C(r) 
4^ 1 



A(r)^ln 
^ ^ 3 \ro 



\-A{r, , 

3 (| + i)(^+i)m. 

0{V~e)- (4.29) 
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Denoting the summand as F{1, e) we see the quantity we 
wish to calculate can be written as 



lim 

e-J-O 



J2Fil,e) + Air)^\ni^ 
li=o ^ ° 



(4.30) 



Eq. (j4.30p is now in a form which is amenable to the 
application of the Riemann sum argument. Using the 
methodology of [11] we expand F{l,e) in the following 
manner: 



F{l,e) = Fo{l)+Fi{xi)^e + ^F{xi,e) 



(4.31) 



where 



Fo(0 = lini^^(/,e) 



L= ((2/ + l)^^^[ln(^o) - V^(U ^o)] 



1 

3 (/ + 1) 



(4.32) 



Here xi = [l + 1/2)-^ and where recall from Eq. 
that Vf) = fcg/8-y/0o. Unfortunately, in this case, it ap- 
pears that we have to calculate the sum of Foil) numer- 
ically. We do this, for a particular space-time, using the 
Levin u transform (see |32| for a detailed discussion of 
this transform). 

Now we turn our attention to Fi(xi). Firstly, we ex- 
pand ^ in the near horizon limit 




2 ( SK^ + r^TQ 
^ 3 I 4K2r3 



(4.33) 



Then we may write F I —= , e | as 



,3/2 



r(i-f^o) 



23/2(^^)1/4 



/3 ^ 4V'o 

X 



1) 



^+1/2 ' 



where 



xf/e + m^rg -|- 1/12 



(4.34) 



(4.35) 



Expanding the last term in Eq. (|4.34p about e = gives 



4-06 



.1/2 



3x/ro\A« 



0(e). 



(4.36) 



We may also expand the second term in Eq. (|4.34p about 
e = to obtain 



xf Ko{axi) 



1 



e3/2 16^K5/4(roa)3/2 



1 



^^^2 (Koixia) {fovox^a^ + 4k (a^ (ro + 2xf ) 

-roxfp)) + 2roXiaKKi{xia) ((4?7iV^ - l)a^ + 4xf/3))} 

+ G(xOe^/^ (4.37) 

where G{xi) is a rather long expression involving combi- 
nations of Bessel functions. 

We now require a similar expansion for the first term, 
which we would expect to cancel with the e — >■ diver- 
gent terms in the above expression, leaving us with a 
finite term that is integrable over x by the Riemann sum 
argument. 

To achieve this we consider the following integral repre- 
sentation for the Whittaker W function: 



2^e 



-2/2 



r(i + ^) 



poo 

Jo 

(4.38) 



We now apply Laplace's method, discussed in Appendix 
IB] for finding an asymptotic solution of this integral as 
1/ oo. Application of this method leads to the following 
expansion for the first term in Eq. (|4.34p : 

xf Ko{axi) 1 

£3/2 16^«;5/4(roa)3/2 

X ^ {xi {Ka{xia) {fUroxfa^ + 4k {a^ (ro + 2xf) 

~roxfl3)) + 2rf)XianKi{xia) ((4mVo - l)a^ +4a;f^))} 
+ G{xiy^ (4.39) 



which we may combine with the expansions (|4.37p and 
(|4.36p . achieving the result 



4?Ao%/e 



3rQK2a:;i 

+\/2xi (rpK + xj) Ki (xia) + roxfy/liKo (xia) 

(4.40) 

One can readily numerically check that AF{xi, e) is both 
0(e) as e — >■ 0, and is a smooth integrable function. 
Therefore, we may apply the Riemann sum argument to 
give the relation 



Sla = Foil) - 



400 A 



Sr^K'^xi 



^3 3/2 -Si/yni" 



+V2xi [r^K + xf) Ki (xia) + roxfy/KKo (xia) 



dx. 
(4.41) 
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Using the identites contained in [3(|, we write down the 
value of the sum Sia in the horizon hmit: 



Sla = Foil) 



a^KTporo (In — 3) 

3%/2 

+ A{r)^ In (^^^+0{e). (4.42) 



Together with Eq. (|4.26p we finally arrive at an expres- 
sion for the sum 5*1 valid in the near horizon limit: 

+ ^(r) — In — H 7 + - ln(e) 

^ ^ 3 \roJ 1920 V "2 



Air) 



2A(A2 + i) 

1 + e" 
°° 2\{\^ + \) 



4^2 

ln(/co(iA)) 



In 



(fco(«A)) +0(e). 



(4.43) 



Stt rQ[g G-gg']analytic = 

34 / 1 X 2^oro(ln(^)-3 

7 - - In {2Krl) ^ — ^ ^ 

1920 V 2 ^ °V 3 



+30fc2p(4/Vo - 27 + 1) + 30 (2fc4o + fc^^) In -i^i^ 
+30(3 - 47)A:4„) - ^ In (ro) + 



(4.47) 



where we recall that fcoo = \Jm?-r\ + 1/12 and Fee' is 
defined in Eq. 

Repeating the procedure for the other derivatives leads to 
the following expressions on the horizon, again splitting 
each term into a numerical and an analytical component 

°° (21 + 1) 

87r^ro[G]„„„ieric = ^ -z (in(i'o) - "0(1 + vq) 



4A 



Renormalized Values 



1 + £2^^^ 

°° 4A 



ln(A;oiA)rfA 



l + e2.AMfco.A)rfA (4.48) 



We are now in a position to substitute the expression 
(14:431) in the definition of {g^^'Cjes'} Eq. (jlje]) and ex- 
pand in the near horizon limit. This procedure yields 



1 



28807r2r4 



0(1). 



{ 1 + 240V'o - 15m2r2(2 + 3TO2r2))} In (e) 

(4.44) 



Comparing this expression with the renormalization sub- 
traction terms p.39p and using the definition of V'o i it is 
straightforward to show that these divergences will can- 
cel exactly with those contained in the subtraction terms, 
giving us the following renormalized value on the horizon: 



I- J analytic 



24 



12i-2 Iv, I ^00 



-12fc^o-ln ( — ) -27-hln(2) 



247^00 



(4.49) 



numeric 



S7:'4[g'* G,u' 

Y^{^l±}1^2klm + u,)^Hu,)) 



8 
1 



3 [l + l) 



'■4)1 + %4l3ii] 
(4.50) 



[5 G;Be']ren = G-ee'] 
where 



I \r,se r 1 

numeric ~r [y ^ -^06' \analytic 

(4.45) 



\numeric 



1=0 



4V'0 1 
3 (l + l) 



-f(2/ + l) 



/(/ + i; 



[ln(i/o) - V(5 + f^o) + 2A 



feoo2A(A2 + i) 

ln(fco(iA)) 



1 



„27rA 



^ e- 

- 2A(A2 + i) 



1 + e 



2n\ 



ln(fco(iA)) 



(4.46) 



4V.i(ln(2roAc)-H) 4V>i 

3 3 " 

rl (5/o"^o - 4ro^(4/^' + 3/^"ro) - 16/«2 
+24kloroiforo + 2k) + (54 - 727)fc^p 
-36/c^ (21n(fcio)-ln(2rgK)) 



1 

288 



40fc?n 



7 



5760 



(-31n(2r2K) -f (67-3)) 



(4.51) 
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•.tt\numeric 

24 



2Aln(fco(»A)) / ro(/^Vo+4/^) 
l + e2^^ V 6 
° 2Aln(fco(*A)) ^ro(/^o+4/^) 

^00 



fco(jA)2 j dA 
A:o(«A)2j dX 

i ^(2Z + + 1) + (m^ + ^i?)r-2]/3oT (4-52) 



1 + e 



2ttX 



6 



i=0 



2V^o 



In(ro) 



+ ^ K (/'"''o - 4ro.^(2/" + 31/"Vo) + le^:^) 



+20A;2o(3fc2o-ro(/"ro + 8K))ln 



2r^K 



-20A;2oro((27 - 5)/"ro + 8(1 + 27)«) + 30(47 - 3)fc^^,] 
288 6 



In {2rl^) + 2(7 - 2)f"rl 



-16(l + 7)Kro] 

Fit 
' 87r2r4 ■ 



40fco'o 7(i_27 + ln(2r^^))r„ 



1920 



(4.53) 



Note the coefficient of (3^ in [5**G;t 



■Mlnumenc COmCS from 

near horizon series expansion of dpoi/dr. It is also found 
that 



[5** G-^tt']ren — [ff'"'" G;rr']r 



(4.54) 



up to terms proportional to R'{r). Of course, due to 
spherical symmetry we must have that 



<f)<f)' \ren 



(4.55) 



Finally we turn our attention to the calculation of 
[g^^G-rr]rem whosc coustructiou follows immediately 
from the above expression. To see this we exploit the 
fact that [M^(x,x']) (the coincidence limit of the regu- 
lar part of the Hadamard expansion for G e ) satisfies the 
inhomogeneous wave equation [l^ | 



{n - - £,R)[W] =-6ui, 



(4.56) 



where for a static spherically symmetric Ricci-constant 
space-time Vi is of the form 



Vl = J^RabcdR'^ 



bed 



1 p r/ab I 
72Q ^abJ^ ~r „ 



1 



1440r4 



{rY"-8/(r/' + l) + /' (8r-4rY') 



+Af + 180toV + 4} . (4.57) 

Now since [g'^^GE-rAren = [g^^W-rr] by definition, we 
have the result that 

W^G-rr]ren — —[g^^G-^tt]ren ~ [g^^ G ■ge]ren 

-[g'^'^G-^Un - {rr? + ei?)[G],e„ - 6i)i (4.58) 

Inserting the expressions calulated for each term on the 
right hand side we arrive at the result that, on the horizon 



\g" G-^rr\ren — [<7*'G;tj]re 



(4.59) 



again up to terms proportional to R'{r). Finally we 
remark that equating the expression obtained here for 
[G]ren with that obtained in [l^ allows one to derive the 
following identity: 



^00 



4A 



o27rA 



ln{koix)dX 



4A 



^00 



„27rA 



ln{koix)dX 



= koo (In(fcoo) - I) - {x, \ + ifcoo) 



x=-l 



ikoa) 



T{x,^+ ifcoo) 



+ ikoo In , — J — ^ , . 



(4.60) 



D. Renormalized Stress Tensor Components 

We begin this section by considering the behaviour of 

{f/)ren-{ft')ren, (4.61) 

near the horizons of a static spherically symmetric black 
hole space-time. This expression, by Eq (|3.6p and the 
fact that Rt* = i?/ (hence TW/ = A^/ by Eq. (E^H)), 
reduces to 

( r }ren }ren — 2^([^? G ■^rr]ren [g G-tt\ren) 



+ (5 ~ OiW^ G-rr']ren — [ff** G-^tt'\ren) 



(4.62) 



As we have already shown that, on the horizon 

\Tr\Ten — and {g'^'G ;ft]ren 5 

we may therefore conclude that 



(4.63) 



on the horizon, r = rg of a general spherically symmetric 
space. 
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Now it is known, that in order for {T^)ren to be fi- 
nite in a freely-faUing frame on the past and future event 
horizons, the fohowing conditions must be satisfied [s^: 

1) \{ft')ren + {f/)ren\ <00. 



2) \{Tg^)ren 

3) 



< oo. 



f 



< OO. 



(4.64) 



Its is straightforward to see that the first two of these 
conditions are satisfied. While Eq. (14.63^ suggests that 
the third condition may be satisfied, the analysis of this 
paper does not allow us to prove this analytically, as we 
cannot rule out terms of the form 0(r — ro)ln(r — tq) 
in the numerator. In order to complete this proof we 
require knowledge of the behaviour of the derivatives of 
(T(*)j.e„ and {T^^)ren as the horizon is approached. If 
they remain finite, then we may say that these quantities 
posses a Taylor series, at least to the first order, about the 
horizon, and hence we may then conclude that condition 
3 holds. 

In Paper II, we calculate both (T^'')ren and (T't*)ren 
on the exterior region, excluding the immediate vicinity 
of the horizons, of a lukewarm black hole. Combining 
these numerical results with the horizon values obtained 
in this paper will allow us to numerically calculate these 
derivatives and thereby draw conclusions on the regu- 
larity of the equivalent of the Hartle-Hawking state for 
lukewarm black holes. We note here that the satisfaction 
of condition 3 can be shown to be equivalent to requiring 
that {Tg^)ren posscss a Taylor series about the horizon, 
by following the method of Morgan et al. [sj] . This will 
also be investigated in Paper II. 

We may now write down expressions for the diago- 
nal components of the stress tensor on a horizon of a 
spherically symmetric black hole space-time. These ex- 
pressions are however quite unwieldy, fortunately with 
the approach adapted here we may easily unite them in 
terms of their sub-components. 



+ |^(i?/ — ^R) — [G]ren 



2vi 
' 8^ 



-4 + X/ (4.65) 



ren — {^(^ 

rr' "1 
^rr' xrr- 



r'\ren \q G'x'r\re'i\) 



87r2 ^ 



(4.66) 



These expressions have the interesting consequence that 
for the minimally coupled case, one need only calculate 
[g^^ Gee'jren, [g'"'' Grr']ren and [G]ren in Order to obtain 
the diagonal components of the stress tensor on a hori- 
zon, with the need for [G]ren disappearing in the massless 
case. 



V. LUKEWARM BLACK HOLES 

We now specialize these results to the case of a luke- 
warm black hole. Lukewarm black holes are a special 
class of Reissner-Nordstrom-de Sitter space-times with 
(Euclidean) line element given by Eq. (|2.7p with metric 
function 



/(r) = 1 - — + — - — 



2M 
r 



(5.1) 



where M, Q are the mass and charge of the black hole re- 
spectively, and A is the (positive) cosmological constant, 
with Q M. For 4M < ^3/A we have three distinct 
horizons, a black hole event horizon at r = r/i, an inner 
Cauchy horizon at r = r_ , and a cosmological horizon at 
r = r^, where 



(5.2a) 
(5.2b) 
(5.2c) 




- i VsM (^1 + ^1- 41/7x73) . 



The fourth root of / is negative and hence nonphysical. 
While the event horizon is formed by the gravitational 
potential of the black hole, the cosmological horizon is 
formed as a result of the expansion of the universe due 
to the cosmological constant 35]. An observer located 
between the two horizons is causally isolated from the 
region within the event horizon, as well as from the re- 
gion outside the cosmological horizon. 
If, as the evidence seems to suggest, the universe pos- 
sesses a cosmological constant [3|| , then it is more natural 
to consider a black hole configuration which is asymptot- 
ically de Sitter than one which sits in an asymptotically 
flat universe. Also given that de Sitter space, in its nat- 
ural vacuum, is awash with radiation to a static observer 
[HI, it seems rather natural that a black hole in a de 
Sitter background would be most comfortable in a final 
configuration in which its event horizon is at the same 
temperature as the surrounding bath. This state of af- 
fairs is realized in the lukewarm case and so the study of 
such a black hole configuration is well motived. In fact 
the lukewarm case has attracted interest recently, as ev- 
idenced in [il,[il[33. 

We shall confine our attention to a single exterior region 
r € [?'/i,fc] which has a regular Euclidean section with 
topology 52 X S'2 [Hi 

In Tables U and |TT] we list the values of the pressure 
{Tg^)ren and the energy density —{T^^)ren for a range 
of valus of the mass of the field on both the event hori- 
zon and cosmological horizons respectively. We consider 
the case of a conformally coupled field of mass m with 
M = O.IL = Q with L = i/Vt. 

As can be seen from the Tables, the energy density on 
the event horizon remains negative as the mass of the 
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field increases, while on the cosmological horizon it re- 
mains positive. The pressure remains positive on the 
event horizon, while its is mostly negative on the cos- 
mogical horizon, as the field mass increases. 
To interpret this result we consider the semi-classical 



m 


Event 


Cosmological 





0.3166133 


-0.0001567 


\L 


0.31458112 


0.0016206 




0.3033364 


-0.0008350 


¥ 


0.2920777 


-0.0013648 


IL 


0.2829577 


-0.0044883 



TABLE I: {Tg^)ren OTi both horizons for various field masses. 



m 


Event 


Cosmological 





-0.0058418 


0.0000407 




-0.0062360 


6.133531 X 10"'' 


hL 


-0.0139295 


0.0007781 




-0.0171130 


0.0019650 


IL 


-0.0117374 


0.0053876 



TABLE IL Energy density on both horizons for various field 
masses. 

field equations: 

If we take the Ag^i, over to the right hand side and con- 
sider it to be a classical stress tensor contribution due to 
the cosmological constant. We see that, since the metric 
has signature (— , and A is positive, the energy 
density of this "stress tensor" and the pressure are both 
negative. One may say that it is this negative pressure 
that leads to infiation. We may conclude that incorporat- 
ing the 'one loop' quantum effects causes an increase in 
energy density on the cosmological horizon while it adds 
to the negative pressure driving infiation for the massless 
conformally coupled case. Similar conclusions for differ- 
ent values of the mass may be drawn from the Tables U 
and mi 



VI. CONCLUSIONS 

The key question in study of quantum field theory on 
black hole spacetimes with multiple horizons, including 
the important class of lukewarm black hole configura- 
tions, is whether it is possible, when the temperature of 
the two horizons is equal, to define an equivalent of the 
Hartle-Hawking state on Schwarzschild space-time which 
is regular on both horizons? In this paper, we have made 
major progress towards answering this question in the af- 
firmative by using the Hadamard renormalization proce- 
dure to calculate expressions for the components of the 
renormalized stress tensor (T^)ren, on the horizons of 
a spherically symmetric black hole spacetime. Making 
use of uniform approximations to the solutions of the 
radial equations, we were able to calculate the value of 
{T^u)ren for a quantum scalar field with arbitrary mass 
and coupling to the geometry and which is in a Hartle- 
Hawking state, on the horizons of a general, spherically 
symmetric black hole spacetime, which is not necessarily 
asymptotically flat. In particular we demonstrated, two 
of the three necessary conditions for regularity: finiteness 

of (ft* + fr'-)ren and of (f/),e„. 

We also made progress in proving the third condition, 
namely, that the \{Tt^ —Tr^')ren\/ (f—Th) be finite by prov- 
ing that {T^'^)ren and {T^*)ren compoucnts are equal on 
the two horizons of our spacetime region. Therefore the 
regularity of the Hartle-Hawking state depends entirely 
on this difference of components possessing a Taylor se- 
ries to the first order about the horizon. To address this 
question, and to find the behaviour of the renormalized 
stress tensor away from the horizon, requires a different 
numerical approach which we will present in Paper II. 
To underline the different nature of these calculation we 
note that to calculate quantities on the horizon it is most 
convenient to separate in a radial direction as only the 
lowest two frequency modes are required there for the 
Hartle-Hawking vacuum. By contrast, away from the 
horizon temporal separation, which is not posssible on 
the horizon since the separation becomes null, is most 
convenient, all modes contribute and we must resort to a 
numerical calculation. In Paper II, we will complete the 
our 'proof of finiteness by providing compelling numeri- 
cal evidence that {Tt*)ren and {Tr^)ren do indeed possess 
finite first derivatives at the horizons. 
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Appendix A: jS}^ and (3^ 

In this Appendix we give expressions for and fi^ 
for the lukewarm black hole configuration. These cal- 
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culations follow along the lines of Appendix B in [l| 
yielding the results, for the event horizon: 



w 



I{r) 



dr + hirc) 



+ 



1 



1 



In (rc - Th) 

H^h - Ti) + hire) 
- b I ln(r/i - r„) 



1 



48(Kr/i2)3/2 
'48y/i^+12k^ { In 



12(27 -l)A;?-r;,(4/,; + /;:) 



(Al) 



which can be shown, numerically, to exhibit a I * be- 
haviour for large I. 

P^i = r ^TriY^(^')dr' + i ln(V.i) + 7 + ln(2ao) 



^ ln(rc - r,,) - ln(r,, - r^) 



1 



ln(r/i - r„) 



which behaves like l^^ for large I. Here 
1 . 1 



(A2) 



& = - {1(1 + 1) + (m^ + ei?)r2 - li(2/; + /^,)) 



+ TT—^ ln(r - r„) - ln(rc - r) 



lir 



a 



{r - ThY 



J 



Finally we note that all the relevant sums over I are 
performed using the Levin u transformation 



Appendix B: Laplace's Method 



In this discussion we draw heavily from Murray |39| . 
Laplace's method is a technique for obtaining an asymp- 
totic expansion of an integral of the form 



lix) 



h(t)x 



(Bl) 



as X — > c». The main concept behind the method is that 
if the function h{t) has its maximum at 0, and ^(O) ^ 0, 
then the dominant contribution to the large x asymptotic 
expansion of I{x) will come from the immediate neigh- 
bourhood of 0. To calculate an asymptotic expansion 
for I{x) as X — >■ cx) it is appropriate to introduce a new 
variable, s, defined by the relation 



h{t) - /i(0) -s^. (B2) 



Hence, the exponential term in Eq. (|B1|) becomes 



^h{t)x _ ^h(0)x^-xs^ 



(B3) 



and so Eq. (jBl[) is of the form 

lix) = / 0(5)6-^"' ds. (B4) 



(A3) By Watson's lemma, an asympto tic p ower series as x 
oo for Eq. (jB4[) may be obtained |39| . 



I{x) 



Following the arguments of jS^ one can obtain an asymptotic expansion of Eq. (|B1[) as x — >■ oo given by 

1 



2tt 



-h"{0)x 



5(0) 



24/i"(0)3x 



1 



125"(0)/z"(0)' + 12/i"'(0)g'(0)/i"(0) + 5(0) (3h'^^\0)h'' (0) - 5/i"'(0) 
144g(4)(o)/i"(0)4 - 210/i"'(0)'/i"(0) (4/1"' (0)5' (0) + 3g{0)h'^^\0) 



-24:h- 



1152/i"(0)6x2 

+21/i"(0)2 (40/?/"(0)V(0) + 5/iW(0) (8/i"'(0)g'(0) + 5(0)/i(4)(0)) + 8giO)h^^Ho)h"' (0) 
"(0)3 (20.g"'(0)/i"'(0) + 15h^^\0)g"{0) + 6h'^^\0)g' (0) + .g(0)/i(5)(0)) + 385g(0)/i"'(0)4] -h 0(x-3)| 



(B5) 
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Higher order coefEcients can be found in 

Now, we shift our focus back to the integral in Eq. (|4.38p . If we introduce a new variable r defined by < !^(r + 1), 
the integral takes the form 



^,(_(,+i)+i„(,+i))^^ + l)-i/2/fo {2^z{t+1)u) dr. 



(B6) 



Setting V — and considering the limit e — >■ allows us to observe, by Eq. (|4.35p . that we require an asymptotic 
expansion of the above integral for v ^ oo. In other words, we may apply the expansion (jBSp to the integral in 
Eq. (jB6[) . Examining this integral we see that, in the notation of Eq. ljBip 



9 



(t) = (r + l)-i/2^o (2v/z(r+l)j/) ; /i(r) = -(r + 1) + ln(T + 1). 



(B7) 



It is clear that /i(t) achieves a maximum value at r = and so we have the following asymptotic expansion for the 
integral in Eq. (|B6l) : 



(B8) 



Jo ' ' M ly I 2Av 

— - \{72vz{2vz - 1) + l)Ko (2y/^) - ASy^il + 2ziy)Ki{2y^)] \ . 



Of course, in this case the function ^(t) is actually a function of both t and v] however we are interested only in 
taking the near horizon limit (z — 0), and in this limit the product vz tends to a constant, so we may treat it as so. 
Next, we consider the integral representation [31] 



(B9) 



Applying Laplace's method to this integral, we obtain the well known large v asymptotic expansion of T[v + i): 

1 1 



T{v + \)7^v''V2^e-'' {l 



2Aiy 1152Z/2 

Combining this expression with Eqs. (|B8p and (|4.38p yields the result 



(BIO) 



1 



1 



24i/ 1552J/2 



I^Ko{2V^) + ^(12j^z - l)Ko{2V^) 



— — — r \(72iyz{2iyz - 1} + l)Xo(2V^) - 48V^(1 + 2zi')Ki(2y^)] 



(Bll) 



Making the replacements 

z^2^oi\ i ^ 2^{a\ + pe^); v = vo- (B12) 

we obtain the expansion (|4.39p . 

It should be noted here that for some parameter sets 
■00 is negative and so i/q is purely imaginary; in this case 
Laplace's method is not applicable, and the method of 
stationary phase must be used instead. In this case we 
consider integrals of the form 




as X ^ oo {x real), and assume that h{t) has a turn- 
ing point at i = c, c € {a,h). The term e"^^*' is now 
purely oscillatory, and as x increases these oscillations 
become more and more rapid, except at the turning point 
c. Near c, the phase xh(t) is approximately constant (or 
stationary) as x — >■ oo. One can expand the integrand 
about this point, in a analogous manner to the analysis 
of Laplace's method, leading to a result similar to Eq. 
(|B5|) to first order. To achieve higher order approxima- 
tions we require the method of steepest descent (see 
for details). Applying this method to Eq. (|B6p when v is 
purely imaginary leads precisely to the result (j4.39p . 
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